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Abstract 

A non-equilibrium extension of Onsager's canonical theory of thermal fluctuations is employed 
to derive a self-consistent theory for the description of the statistical properties of the instanta- 
neous local concentration profile n(r, t) of a colloidal liquid in terms of the coupled time evolution 



equations of its mean value n{r,t) and of the covariance a{r,r';t) = 6n{r,t)6n{r' ,t) of its fluc- 
tuations 5n{r,t) = n{r,t) — n{r,t). These two coarse-grained equations involve a local mobility 
function 6(r, t) which, in its turn, is written in terms of the memory function of the two-time cor- 



relation function C{r,r';t,t') = 6n{r,t)6n(r',t'). For given effective interactions between colloidal 
particles and applied external fields, the resulting self-consistent theory is aimed at describing the 
evolution of a strongly correlated colloidal liquid from an initial state with arbitrary mean and 
covariance rf{r) and a^{r,r') towards its equilibrium state characterized by the equilibrium local 
concentration profile n^'^(r) and equilibrium covariance a^'^{r,r'). 

This theory also provides a general theoretical framework to describe irreversible processes as- 
sociated with dynamic arrest transitions, such as aging, and the effects of spatial heterogeneities. 



PACS numbers: 05.40.-a, 64.70.pv, 64.70. Q- 
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I. INTRODUCTION 



In this paper a non-equilibrium generalization is presented of the self-consistent general- 
ized Langevin equation (SCGLE) theory of colloid dynamics 1, 2!, and of its recent adap- 
tation as a theory of dynamic arrest 4], with the purpose of describing non-equilibrium 
diffusive phenomena in general, and irreversible aging processes associated with the glass and 



the gel transitions 



-|9| in particular. This generalized theory is based on a non-equilibrium 



extension of Onsager's canonical theory of thermal fluctuations. The resulting theory con 



tains, for example, the fundamental equation of dynamic density functional theory 



10| as a 



particular limit, whereas in ot 
early spinodal decomposition 



ler limit one can recognize the basic equation of the theory of 



11| . A practical and concrete use of the resulting general the- 



12| with a quantitative apphcation 



ory of colloid dynamics is illustrated in a related paper 
to the prediction of the aging processes occurring in a suddenly quenched colloidal liquid. 
The dynamic properties of colloidal dispersions has been the subject of sustained interest 



for many years [13 



15| . These properties can be described in terms of the relaxation of 



the fluctuations 6n{r,t) of the local concentration n{r,t) of colloidal particles around its 
bulk equilibrium value n = N/V. The average decay of 6n{r, t) is described by the two- 



time correlation function F{k,T;t) = V ^(5n(k, t -|- r)(5n(— k, t) of the Fourier transform 
6n(]<.,t) of the fluctuations 6n{r,t), whose equal-time limit is S{k;t) = F{k,T = 0;t) = 



V^^6n{]<.,t)6n{—]<.,t). We shall refer to the time r as the correlation time. If some external 
(or internal) constraints that kept a system at a certain macroscopic state are broken at the 
(evolution) time t = the system relaxes spontaneously, searching its new thermodynamic 
equilibrium state. If the end state, however, is a glass or a gel, one refers to t as the waiting 
or aging time The evolution of S{k;t) and F{k,T;t) as a function of the time t 

characterizes the non-equilibrium evolution of the system, and its theoretical understanding 
is a major fundamental challenge. 

If the system is a fluid and it has fully relaxed to its thermodynamic equilibrium state, 
the properties above no longer depend on t, i.e., F{k,T;t) = F{k,T) and S{k;t) = S{k). 
The equilibrium stationary correlation function F{k, r) is then referred to as the intermedi- 
ate scattering function, and its initial value S{k) as the equilibrium static structure factor. 
These properties can be measured by a variety of experimental techniques, including (static 
and/or dynamic) light scattering 6,[3,[3|. S{k), being an equilibrium property, is amenable 



2 



to theoretical calculation using statistical thermodynamic methods The fundamental 
understanding of F{k, r), on the other hand, requires the development of theoretical meth- 
ods to describe the correlations of the local concentration fluctuations, and a number of such 



13 



la, UM- One of them 



approaches have been proposed for their theoretical calculation 
has been developed within the last decade and is referred to as the self-consistent gener- 
alized Langevin equation (SCGLE) theory of colloid dynamics Q, y, Q-Q- This theory 
has been recently applied to the description of dynamic arrest phenomena in several spe- 
cific colloidal systems that include mono-disperse suspensions with hard-sphere interactions. 



moderately soft-sphere and electrostatic repulsions, short-ranged attractive interactions, and 



21 



26 



model mixtures of neutral and charged particles [3|, 

In spite of the long tradition in the study of glasses |27H29[|. until recently the only 
well-established and successful theoretical framework leading to first-principles quantitative 
predictions of the dynamic properties of colloidal liquids near their dynamic arrest transi- 
tion was the conventional mode coupling theory (MOT) of the ideal glass transition 291-132 1. 
Many of the predictions of this theory have been systematically confirmed by their detailed 
comparison with experimental measurements in model colloidal systems 33-|40]. In this 
context, we can mention that the more recently-developed SCGLE theory of dynamic ar- 
rest leads to similar dynamic arrest scenarios as MOT 3|, |2l| for several specific (mostly 
mono-disperse) systems, although for colloidal mixtures differences may appear in some 
circumstances, as reported in Refs. 24-2^. 

An important common feature of both theories in their current status is that they are 
able to predict the regions of the control parameter space in which the system is expected 
to be dynamically arrested, i.e., they predict what we refer to as the "dynamic arrest phase 



diagram" of the system 



25 



26|. While it is important to pursue the application of these 



two theories to specific idealized or experimental model systems and to compare their pre- 
dictions, it is also important to attempt their extension to the description of the detailed 
non-equilibrium processes leading to dynamically arrested end states. Aging effects, for ex- 
ample, should be a fundamental aspect of the experimental and theoretical characterization 
of these non-equilibrium states. These preoccupations have been addressed in the field of 
spin glasses, where a mean-field theory has been developed within the last two decades 41 |. 
The models involved, however, lack a geometric structure and hence cannot describe the 
spatial evolution of real colloidal glass formers. Although experimental studies US and 
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computer simulations 42|-|44j have provided important information about general proper- 



ties of aging, until now no quantitative first-principles theory is available to describe the 
irreversible formation of structural glasses. 



About a decade ago Latz 



45| attempted to extend MCT to describe the irreversible re- 



laxation, including aging processes, of a suddenly quenched glass forming system. A major 
aspect of his work involved the generalization to non-equilibrium conditions of the conven- 
tional equilibrium projection operator approach \^ to derive the corresponding memory 
function equations in which the mode coupling approximations could be introduced. Sim- 
ilarly, De Gregorio et al. 4^ discussed time-translational invariance and the fluctuation- 
dissipation theorem in the context of the description of slow dynamics in system out of 
equilibrium but close to dynamical arrest. They also proposed extensions of approximations 
long known within MCT. Unfortunately, in neither of these two theoretical efforts, quanti- 
tative predictions were presented that could be contrasted with experimental or simulated 
results in specific model systems of structural glass-formers. 

The present work is aimed at extending the SCGLE theory of dynamic arrest to non- 
equilibrium conditions. This paper contains the proposal of such general theory while the 
accompanying paper (paper 11) reports a concrete quantitative application. The general 
theory proposed here consists of the time evolution equations for the mean value and for 
the covariance of the instantaneous local concentration profile n(r, t) of a colloidal liq- 
uid coupled, through a local mobility function b{r,t), with two-time correlation function 



C{r,r';t,t') = 6n{r,t)6n{r' ,t') of the local concentration fluctuations. A set of well-defined 
approximations in the memory function of C(r, r'; t, t') leads to the non-equilibrium extension 
of the self-consistent generalized Langevin equation theory of colloid dynamics to spatially 
non-uniform and temporally non- stationary systems. The resulting theory is applied in II 
to the description of aging effects in a specific model glass forming colloidal liquid. 

In contrast with MCT, the SCGLE theory does not involve the assumption of an under- 
lying Hamiltonian (or any other microscopic) level of description, nor the use of projection 
operator techniques. Instead, it is based on what we refer to as Onsager's canonical theory 
of equilibrium thermal fluctuations. Since the description of thermal fluctuations and relax- 
ation processes can be approached from a bewildering number of theoretical perspectives, 



involving a diversity of issues, approaches, aims, methodologies, and nomenclature 



3-501, 



it is necessary to state that in this work for "Onsager's theory" we mean the general and 
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fundamental laws of linear irreversible thermodynamics and the corresponding stochastic 



;heory of thermal fluctuations, as stated by Onsager 



53 



51 



54l |. respectively, with an adequate extension [55 



52| and by Onsager and Machlup 



56| to allow for the description of 



memory effects. 

Viewed as a theory of fluctuations, Onsager's theory refers to systems in thermodynamic 
equilibrium, and hence, assumes stationary conditions. Thus, generalizing the SCGLE the- 
ory of colloid dynamics to non-equilibrium calls for an extension of Onsager's theory to 
non-stationary non-equilibrium conditions, outside the so called "linear regime", where its 



validity has been universally tested 



48| . Such an extended Onsager's theory is discussed 



in detail elsewhere 57|, and here we only provide a brief summary (see Sec. HIl below). In 
essence, however, this extension consists of the assumption that the t-dependent irreversible 
evolution of a system towards its stable equilibrium state proceeds as a virtually continuous 
sequence of non-equilibrium, but momentarily stationary, states. The main objective of the 
present paper is then to apply this extended canonical theory as a fundamental framework 
in which to discuss the dynamics of a colloidal suspension that evolves irreversibly towards 
its equilibrium state. Such application is the subject of Sec. IIIII 

According to this program, our discussion will involve two distinct levels of generality. 
The first corresponds to the rather abstract and most general description provided by On- 
sager's extended theory in terms of a set of macroscopic state variables, generically denoted 
by (ai, 02, flAf) = a, as reviewed in the following section. The second corresponds to 
the description of diffusive processes in colloidal dispersions, where the abstract objects in 
Onsager's theory take a concrete meaning. Bridging these two levels of discussion requires 
that we identify the specific correspondence between the abstract concepts in Onsager's 
theory and the concrete concepts pertaining to the other more specific level. For example, 
the abstract state variables will be identified with Ni/ AV , the number concentration of 
particles in the ith cell of an imaginary partitioning of the volume occupied by the colloidal 
system in M cells of volume AV . In the continuum limit, the components of the state vector 
a(t) then become the local concentration profile n(r, t) and the fundamental thermodynamic 
relation S = S\a] (which assigns to any point a of the thermodynamic state space a value of 
the entropy [58]) will be identified with the functional dependence of the free energy on the 



ocal concentration profile employed, for examp 



59| or in its more recent dynamic version 
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e, in the classical density functional theory 



60 



61| . For completeness, the structure of 
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this thermodynamic framework is reviewed in an Appendix. Finally, in the last section we 
summarize the main conclusions of the present work. 



II. GENERALIZED ONSAGER THEORY 



In this section we summarize the main features of the extension of Onsager's theory 



to non- stationary non-equilibrium states presented in detail in Ref. [57|]. Thus, consider 
a system whose macroscopic state is described in terms of a set of M extensive variables 
aj(t), i = 1, 2, M, which we group as the components of a M-component (column) vector 
a(t). The fundamental postulate of this generalized theory is that the dynamics of the state 
vector a(t) may be represented by a multivariate stochastic process which is globally non- 
stationary, but that within any small interval of the evolution time t may be regarded as 
approximately stationary. This local stationarity approximation is then complemented with 
the assumption that the mean value a(t) is the solution of a generally nonlinear equation, 
represented by 

^ = 7^[a(t)], (2.1) 
whose linear version in the deviations Aa(t) = a(t) — a^'^ from an equilibrium value a^'' reads 

^^ = -£[a^^]-^[a^«]-Aa(t), (2.2) 

with C and S being MxM matrices and with the symbol "■" indicating the corresponding 
matrix product. The matrix £[a'^'^] is referred to as the kinetic matrix, related with the 
vector of "fluxes" 7^ [a^«] of Eq. (El]) by Cla^"] = - {dfR [a] /9a) ■ E'^ [a'^"]. 
On the other hand, E [a] is the thermodynamic matrix, defined as 

in which the function 5 [a] determines the dependence of the entropy on the components of 
the vector a, i.e., S* = ^[a] is the fundamental thermodynamic relation of the system, and 
hence, Fj\s\ = k~^^ {dS[a]/daj) is the conjugate intensive variable associated with aj. One 
should notice that Eq. can be written as dAa{t)/dt = £[a^5] ■ AF(t), where AF(t) = 

F[a(t)] — F^'' is the macroscopic deviation of the vector F[a(t)] of intensive parameters from 
its equilibrium value F*^^ = F[a^'']. This relaxation equation is immediately recognized as 
the classical format of the linear laws of irreversible thermodynamics. 



From these premises a time-evolution equation for the MxM covariance matrix 



6a{t)6a>{t) can be derived [57|, which reads 
da{t) 



dt 



£[a(t)] ■ £ [a(t)] ■ a{t) - a{t) ■ £ [a(t)] ■ C^si{t)] + {C[si{t)] + /:^[a(t)]) . (2.4) 



This equation may be regarded as a simple extension of the equation of motion for the 
covariance involved in the conventional Onsager theory (see, for example, Eq. (1.8.9) of 
Ref. |48|), in which the matrices /^[a"^^] and £^ [a'^''] are replaced by >C[a(t)] and £^ [a(t)]. 
The detailed arguments to see that this is the proper manner to extend Onsager's result to 
non-stationary processes can be found in Ref. [57|. 

Thus, if two essential pieces of information were available, namely, the fundamental ther- 
modynamic relation S = S[a] and the state-dependence of 71 [a], then Eqs. (12. ip and (12. 4p 
would constitute a closed system of equations for the mean value a(t) and the covariance 
cr{t). These are essentially the first and second moments of the 1-time probability distribu- 
tion Pi(ai,ti) that the state vector a has the value ai at the time t = ti. Let us stress 
that the full knowledge of this probability distribution is equivalent to the full knowledge 
of the macroscopic state of the system, even though for many purposes one may only be 
interested in one or some of its moments. For example, the probability distribution of the 
thermodynamic equilibrium state is fully determined by the Boltzmann-Planck postulate to 
be 48l. l62l| Pi fa. t) = P^''[a] = exp [{S[a] — S[a.'^'^]) /ks], whose mean value and covariance 



are 



58 



63| a(t) = a^'^ and a{t) = cr^'^ = £ ^ [a^'^]. One may stretch this concept, and 



introduce the assumption that the non-equilibrium evolution of the system can be described 
approximately by Pi(a, t) = P^'"'^'\a,t) = exp(5'[a] — S\a(t)])/kB, whose covariance is given 
by a(t) = a^^'^ \t) = £^^ [a(t)]. We shall refer to this as the local equilibrium approxima- 
tion, and an idealized time- dependent process that satisfies this approximation at any time t 
shall be referred to as a quasi-static process. We must emphasize that the present extended 
Onsager's theory is, of course, NOT based on this approximation. 

If the goal were to fully determine Pi(a, t), in principle, one could attempt to write 
the time-evolution equations corresponding to the higher-order moments, thus constructing 
an infinite hierarchy of equations for all such moments. Alternatively, also in principle, 
one could attempt to write the time-evolution equation for Pi(a, t), from which one could 
determine the time-evolution of all the moments. Our intention, however, is not to follow any 
of these strategies, nor to assume that the stochastic process is Gaussian, so that the first two 
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moments above will suffice to fully determine Pi (a, t). In fact, we are not actually interested 
in determining Pi (a, t) at all. Instead, our aim is to use the two general equations above for 
a(t) and a{t) as the fundamental framework in which to introduce approximations that lead 
to a closed system of equations for these directly measurable properties, at least in specific 
and concrete cases, as in the colloidal context discussed below. For this purpose, rather 
than analyzing the higher-order moments of Pi(a, t), we consider the two-time correlation 
function, i.e., the second moment of the two-time probability distribution P2(ai, ti; a2, ^2), 
as one aspect of the properties of the thermal fluctuations (5a(t + r) = a(t + r) — a(t) around 
the non-stationary mean value a(t), within the local stationarity approximation. 

Thus, the second fundamental postulate of the generalized Onsager theory is that the 
locally stationary fluctuations 6a{t + r) can be described by a mathematical model that we 



refer to as a generalized Ornstein-Uhlenbeck process, discussed in Ref. [55|, which in the 
present context is defined by the most general linear stochastic differential equation with 
additive noise, which has the following general structure 

'^^^'^1 + ^'^ = -c^[a(t)]-[a(t)]-i-5a(t + r)- r dr'^[T -T'-Mt)]-[a{t)r' ■ Sait + r') + i{t + T), 
dr Jo 

(2.5) 

in which the stochastic vector f(t + r) is assumed stationary but not necessarily Gaussian or 
5-correlated, the matrix u[a] is antisymmetric, a; [a] = — a;^[a], and the memory matrix 7[r; a] 
satisfies the fluctuation-dissipation relation 7[r; a(t)] = 7^[— r; a(t)] = < f(t + r)f ^(t + 0) >. 
From this generalized Langevin equation one then derives the time-evolution equation for 



the non-stationary time-correlation function C(r; t) = 6a{t + r)5a'l"(t), which reads 

= -uim] ■ ^~\t) . C(r; t) - j^^ dr'7[r - r'; a(t)] ■ a-\t) ■ C{t'- t), (2.6) 

and whose initial condition is C(r = 0;t) = cr{t). This equation describes the decay of 
the correlation function C{T\t) with the "microscopic" correlation time r, after the system 
has evolved during a "macroscopic" evolution time t from an initial state described by 
a*^ = a(t = 0) and cr° = o"(t = 0), to the "current" state described by a(t) and cr(t). 

Notice that this equation involves a{t) explicitly, and a(t) implicitly through the matrices 
c<;[a(t)] and 7[r; a(t)]. Thus, besides requiring the actual values of a(t) and cr()f:), this equation 
also requires the information on the state dependence of the matrices a; [a] and 7[r; a]. This 
information, however, must be closely related with the kinetic matrix £[a]. To establish 
such relationship, notice that if the system has reached a thermodjTiamic equilibrium state, 
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in which Pi(a,t) = P^''[a], a(t) = s!"^, and a{t) = cr'^" = S''^ [a^"], then Eq. ([23]) may be 
written, in the so-called "Markov" limit, as 

= -£[a^^] ■ S [a^^] ■ 5a(r) + f(r), (2.7) 

with C[a] defined as 

poo 

C[a] =cu[a] + / rfr7[r;a]. (2.8) 
Jo 

Eq. (12. 7p . however, is the linear stochastic equation with additive white noise of the Onsager- 



Machlup theory of equilibrium fluctuations 53|, |5J|. According to Onsager's regression 



hypothesis, this equation must be identical, except for the additive white noise f(r), to the 
phenomenological relaxation equation in Eq. (12. 2p . This requires that the deflnition of the 
matrix C[a] in Eq. (12. 8 p above must be consistent with the phenomenological definition 
£[a^«] = - (dn [a] /da)^^^,, ■ [a^i] beneath Eq. (Q- 

The results above then state that the kinetic matrix C[a] may be obtained either by lin- 
earizing the non-linear phenomenological relaxation equation (12. ip if this equation is known 
a priori, or by means of the general relationship in Eq. (12. 8 p if the matrices a;[a(t)] and 
7[r;a(t)] can be determined by independent arguments, as we propose here in the context 
of colloid dynamics. In general, the antisymmetric matrix a;[a(t)] represents conservative 
(mechanical, geometrical, or streaming) terms, and its determination in specific contexts is 
relatively straightforward. In contrast, the memory matrix 7[r;a(t)] summarizes the effects 
of all the complex dissipative irreversible processes taking place in the system. Its exact de- 
termination is perhaps impossible except in specific cases or limits; otherwise one must resort 
to approximations. These may have the form of a closure relation expressing 7[r;a(t)] in 
terms of the two-time correlation matrix C(r; t) itself, giving rise to a self-consistent system 
of equations, as we illustrate in the application that follows. 

As a final observation, let us mention that throughout the previous discussion we have 
assumed that the variables aj(t) represent extensive state variables. In reality, one could 
also describe the state of the system in terms of any combination of extensive and inten- 
sive variables. The choice depends on the convenience, given the macroscopic conditions 
imposed on the system. Using only extensive variables, for example, is the most convenient 
choice if the system is subject to isolation conditions. If, however, the system is in contact 
with a thermal reservoir, the temperature, rather than the internal energy, may be a more 
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convenient variable. On the other hand, if the external constraints (isolation, contact with 
thermal reservoirs, applied external fields, etc.) are time-independent, the time-evolution 
equations for the mean value a(t), for the covariance a{t) and for the correlation function 
C{T]t) (Eqs. (12.11) . (12. 4p . and (12.61) ) will describe the spontaneous relaxation of the system 
toward the corresponding equilibrium state. We may, however, also consider the possibility 
that these constraints vary in time in a programmed manner. In this case, the time-evolution 
of the parameters describing these constraints (for example, the overall density of the system 
or the temperature of the heat reservoir) may be prescribed, rather than determined as the 
solution of any time-evolution equation such as Eq. (12. ip . and the time-evolution equations 
for the covariance and the correlation function (Eqs. (12. 4p and (12. 6p ) will describe the non- 
equilibrium response of the system to these forced time- dependent macroscopic constraints. 

III. APPLICATION TO COLLOID DYNAMICS 

In this section we discuss the general problem of the diffusive relaxation of the local 
concentration of colloidal particles in the absence of hydrodynamic interactions but which 
interact through pairwise direct forces represented by the effective pair potential n(r,r'). 
Thus, let us consider a dispersion of such colloidal particles of mass m in a volume V 
which, in the absence of external fields, has a uniform bulk number concentration ub = N/V . 
In the presence of a conservative static external field that exerts a force F'^^*(r) = — V'?/'(r) on 
one particle located at position r, the mean local concentration profile of colloidal particles, 
rr(r, t), will evolve in time from some initial condition n(r, t = 0) = n^{T), towards its stable 
thermodynamic equilibrium value n^'^(r), while the covariance (j(r, r'; t) = 6n{r, t)6n{r', t) of 
the fiuctuations 5n(r, t) = n(r, t) — n(r, t) will evolve from an initial value cr°(r, r') to a final 
equilibrium value cr^'(r, r'). The initial values n°(r) and o"''(r,r') are, of course, arbitrary, 
whereas the final equilibrium mean and covariance, n'^''(r) and cr^^(r, r'), are univocally dic- 
tated by the external constraints imposed on the system (isolation, contact with reservoirs, 
etc.) and by the external field ip{v), according to the second law of thermodynamics. We 
open this section with a brief reference to the specific thermodynamic framework in which 
this problem is embedded. Given its conceptual importance, we provide additional details 
on this topic in Appendix A. In the rest of this section we elaborate the dynamic aspects as 
a concrete application of the generalized Onsager theory just reviewed. 
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A. Thermodynamics of fluids in spatially inhomogeneous states. 



The most fundamental thermodynamic ingredient in the apphcation of this general the- 
ory is the fundamental thermodynamic relation (FTR) S = S[a], which assigns a value of 



the entropy at any possible values of the state variables (ai, 02, ) = a SSj. In practice, 



however, we only need the first and second derivatives of S'[a], which define the intensive 



'3^ 



parameters -Fj[a] = k]^ {dS[a]/daj) and the thermodynamic matrix <6^y[a] = —{dFi[a]/da 
To specify the variables (ai, 02, oa/) = a of our problem, let us first mentally partition the 
volume ^ in a number C of smaller portions (or cells), whose internal energy, particle num- 
ber, and volume, we denote by E^'^\ N^'^'' and V^^\ respectively, with r = 1,2,...,C. Then, 
the fundamental thermodynamic relation of this system reads S = 5'[E,N, V], where E, N 
and V are C-dimensional vectors with components E^'^\ N^'^\ and V^"^^ (r = 1, 2, C). For 
the sake of simplicity let us assume that the volumes V'-''^ are all equal, V'-''^ = AV = V/C, 
and remain fixed, so that only the variables [E,N] are needed to define a thermodynamic 
state. This corresponds to the selection a = [E,N]. 

Just like in ordinary classical thermodynamics, under some circumstances one may prefer 
to express the FTR not in terms of the variables [E,N], which involve the local internal 
energy E^'^\ but in terms of the particle number profile N and of some form of "local 
temperature". Such representation is most convenient under conditions in which the 
particle system is in contact with a thermal reservoir (in our case the supporting solvent) that 
keeps temperature constant and uniform. Under these circumstances the chemical equation 
of state can be expresses as the dependence [(3^'] N] of the local electrochemical 

potential on the profile N and of the thermal reservoir parameter as it is explained in the 
appendix. In the general expression for N] in Eq. (13. ip below the explicit reference 

to the parameter is omitted. 

Let us notice that the discussion above is independent of the spatial resolution employed 
to describe the non- uniformity of the distribution of matter and energy, i.e., on the number C 
of cells in which we mentally partitioned the total volume V. Since this is a mere informatic 
concept (the cells are not meant to represent macroscopic subsystems), one can take the limit 
of maximum resolution AV — )■ (or C — > 00). Although no new concepts arise in taking 
this limit, the notation and the nomenclature change somewhat. For this, let us define the 
vectors n and e whose components are the local particle number density n*-'"-' = N^"^^ / AV 
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and the local energy density e^"^^ = E^^^AV, and whose average remain finite in this limit. 
Second, rather than labeling the cells with the discrete index r, running from 1 to C, now 
we label them with the position vector r of their centers. In the limit of vanishingly small 
cells, r varies continuously in the volume V, and hence, the vector components n^''''^ and e^''^ 
with r = 1, 2, C now become the functions n{r) and e(r) of the position vector r G V. As 
a consequence, what used to be an ordinary function of the vectors N and E, such as the 
entropy, now becomes what is called a functional of the functions n{r) and e(r). For example, 
the local electrochemical potential [f3^; n] of the particles at cell r now becomes an 

ordinary function of the position vector r, and a functional of n{r). This dependence will be 
indicated as /i[r; P^, n] or simply as /i[r; n]. Of course, the ordinary derivative of a function, 
such as the thermodynamic matrix S^'''''^''^[N; /3] = ((9/3/i*^''^[N; /3]/(9A^'^'''^) in Eq. f|T5|) of the 
appendix, now becomes the functional derivative (5/3/i[r; n]/6n{r')). In the continuum limit, 
we must also replace {AV) by a volume integral J (Pr on the vector r, and 5^.^ / AV by 
the Dirac delta function, 5(r — r'). Clearly, "matrices" such as u^"^'^'^ now become functions 
of the two position vectors r and r', and matrix products now become convolutions. 

With this notation, let us now write the most general expression for the local electro- 
chemical potential yu[r; at position r in units of the thermal energy fc^T = namely 

/3/i[r;n] = /3/i'"[r; n] + /3^(r) (3.1) 
= + Inn(r) - c[r; n] + /3^(r). 

In this equation ip{v) is the potential of the external field acting on a particle at position 
r. The first two terms of this definition of yU*"[r;r;,], (/3yU*(/3) + Inn(r)), are the ideal gas 
contribution to the chemical potential, whereas the term — c[r; n] contains the deviations 
from ideal behavior due to interparticle interactions. 

Using Eq. (13. ip . the thermodynamic matrix £^[r, r';?7,] = [5/3yu[r; n]/(5n(r')] can then be 
written in general as 

£[r, r'; n] = S{r - r')/n{r) - c^^) [r, r'; n], (3.2) 

with c*^^^ [r, r'; n] = {5c[r;n\/Sn{r')) being the functional derivative of c[r; n] with respect 
to n(r'), referred to as the direct correlation function. On the other hand, the covariance 
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matrix cr(r, r') = Sn{r, 0)Sn{r', 0) can be written in terms of the total correlation function 
/i*^^)(r, r') as 



a{r, r') = n{r)6{r — r') + n{r)n{r')h^'^\r, r'). (3-3) 

The matrices S[r, r'; n] and cr(r, r') are not in general related to each other. It is only when 
they are evaluated at the equilibrium local concentration profile n^''(r) that they are related 
to each other by the second equilibrium condition in Eq. (fTSjl of the appendix, which in the 
present notation reads 

j (irV"«(r, r')£[r', r"; n'^] = 5{v - r"). (3.4) 

Using Eqs. (13. 2p and (13. 3p . one can immediately see that this equation is equivalent to the 
well-known Ornstein-Zernike equation 16| 



h{r, r') = c(r, r') + J rfV"c(r, r"X^{r")h{r" , r'), (3.5) 

where c(r, r') and /i(r, r') are, respectively, the equilibrium value of c*^^^(r, r') and /i*^^^(r, r'). 

In conclusion, the thermodynamic matrix S[r,r']n{t)] evaluated at an arbitrary state 
n(t) is fully determined by the chemical equation of state. Its equilibrium value, £[r, r'; n^'^], 
determines the covariance a'^'^{r,r') of the equilibrium distribution by means of Eq. (13.41) . 
which is equivalent to the Ornstein-Zernike equation above. The time- dependent covari- 
ance cr(r, r'; t) of an arbitrary non-equilibrium state, however, cannot be determined in this 
manner, unless the local equilibrium approximation (i.e., the quasi-static limit) is assumed 
to be valid. Thus, in general we need an independent, non-thermodynamic condition, to 
determine this important property, and this is the main subject of the following subsection. 



B. (Irreversible) Time evolution of n(r,t) and cr(r,r';t). 



Let us start by writing the analog of Eq. (12.11) . The macroscopic diffusive relaxation of 
the local concentration n{r, t) of colloidal particles is described by the most general non- 
linear but spatially and temporally local diffusion equation provided by Fick's law, which 



reads 



48 



49| 



dn{r, t) 
dt 



■ 6(r, t)n{r, t) V/3/i[r; n{t)]. 



(3.6) 
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In this equation Dq is the diffusion coefficient of the colloidal particles in the absence of 
interactions between them and b{r, t) is a local reduced mobility, to be specified later, which 
describes the frictional effects of the direct (i.e., conservative) interactions between particles, 
as deviations from the value b{r,t) = 1. 

We may now linearize this equation around the equilibrium profile n'^''(r), to get the 
analog of Eq. (12. 2p . from which we can identify the "matrix" C[r, r';n{t)] of Onsager kinetic 
coefficients as 



- £[r, r'; n{t)] = D°V ■ n(r, t) 6(r, t)V5(r - r'). (3.7) 
Using Eq. (12. 4p we can then write the relaxation equation for a{r, r'; t) as 

^^^^^ =D°V ■ rl(r, t) 6(r,t)V J dr,£[r,r,;n{t)Hr2,r'-t) 

+ D^V ■ n{r', t) b{r', t)V j dv^Sir', r^; n(t)]a(r2, r; t) (3-8) 

- 2D°V • n(r, t) 6(r, t)V6{r - r'). 

Let us now describe the fluctuations 6n{r,t + r) = n{r,t + r) — n{r,t) of the local 
concentration at position r and time t + r around the mean value n(r, t) within a microscopic 
temporal resolution described by the time r. The assumption of local stationarity means 
that in the time-scale of r, n{r,t) is to be treated as a constant. We may add the spatial 
counterpart of this simplifying assumption. Thus, we write the fluctuations as 5n(r + x, t + 
r) = n(r+x, t+r)— rr(r, t), where the argument r of n(r, t) refer to the macroscopic resolution 
of the measured variations of the local equilibrium profile, whereas the position vector x 
adds the possibility of microscopic resolution in the description of the thermal fluctuations. 
Defining the fluctuations as the deviations of the microscopic local concentration profile 
n{r + X, t + r) from the mean value n{r,t) indicates that, within the microscopic spatial 
variations described by the position vector x, n(r, t) must be treated as a constant. To a 
large extent, this is equivalent to recover the partitioning of the system in cells of a small but 
finite volume AV, and assume that variations from cell to cell are described by the vector r, 
whereas variations within cells are described by the vector x, and that within the intra-cell 
scale, the system can be regarded as uniform and isotropic. Under these conditions, the 
covariance a{r + x, r + x'; t) may be written as cr(| x — x' |; r, t), and in terms of its Fourier 
transform cr(A;; r, t), as 
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a(| x-x' \;r,t) = ^ J d'ke-''^<---'^a{k;r,t). (3.9) 
In this manner, Eq. (13. 8 P may be re-written as 



'^'^^^^'''^^ = - 2eD^n{r, t)b{r, t)S{k; n{r, t))cx{k; r, t) 
+ 2k'^D^n{r,t) 6(r,t), 

where S{k;n{r,t)) is the FT of S{\ x — x' \;n{r,t)), defined as the thermodynamic matrix 
evaluated at a uniform concentration profile with a constant value given by the local and 
instantaneous concentration rl(r, t) at position r and time t. 

C. Relaxation equation for C{t,t') 

The description of the fiuctuations 5n(r + x, t + r) = n(r + x, t + r) — n(r, t) with the 
temporal and spatial resolution described by the time r and position vector x cannot be 
obtained by simply linearizing the macroscopic version of Pick's diffusion equation above. 
Instead, one has to consider a generalized version of Pick's law, which contains Eq. (13. 6p as 
its macroscopic limit. Such an extension reads 



dn{r^ t) 
dt 



D°V • ^ dt' j d^r' b[r -r']t- t']n{r' , t')VPfi[r'; n{t% (3.11) 



where 6[r; t] is a time-dependent local mobility that must yet be specified. 

This generalized diffusion equation may be derived rather simply by complementing the 
continuity equation, 

^ = -V.j(M). (3.12) 

with a constitutive relation constructed at the level of the particle current. We require that 
the friction force on the particles in the neighborhood of position r must be equilibrated 
by the osmotic force — V/i*"'[r; n] and by the external force — V\E'(r) on each particle, both 
included in — V/i[r;n], so that 

C°j(r,t) + ^ dt' j d\' AC[r-r'-t-t']-i{r,t) = -n{r,t)Vfi[r;n{t)], (3.13) 

The friction force per unit volume on the left hand side of this equation is the sum of 
the friction due to the supporting solvent, (^}{r,t), and the frictional effects due to the 
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interactions between the colloidal particles themselves, (AC)j(r,t). The latter, however, is 
assumed to be in general spatially and temporally nonlocal. The solution of this equation 
for j(r, t) can be written as 

j(r, t) = -D^ dt' J rfV6[r - r; t - t']n{r', t')V/3fi[r'; n{t')], (3.14) 

where is the free diffusion coefficient, defined here as = ksT/C^^, and where the 
spatially and temporally non-local mobility kernel 6[r — r'; t] is defined in terms of the 
memory function AC*[r — r'; t — t'] = AC[r — r'; t — t']/C° as the solution of the equation 



b[r - r'; t] = 6{r - r')26{t) - dt' j dV"AC*[r - r"; t - t']6[r" - r'; t']. (3.15) 

Using Eq. f l3.14p in the continuity equation fl3.13p finally leads us to Eq. (13. lip , which 
reduces to Eq. (13. 6 p when the generalized mobility kernel 6[r — r'; t — i:'] is approximated by 
its spatially and temporally local limit. 



h[v -r';t- t'] = b{r, t)5{r - r')26{t - t'), (3.16) 



where 



//•CXD 
die J dr 6[x,r;r,it] (3.17) 

with 6[x, r; r, t] = b[{r + x) - r; (t + r) - t] . 

We can now proceed to identify the elements of Eq. (12. 5p corresponding to our prob- 
lem. In the present case, the corresponding antisymmetric matrix w[a(t)] vanishes due to 
time-reversal symmetry arguments [55]. We can then write the matrix 7[r;a(t)] as the non- 
markovian and spatially non-local Onsager matrix implied by the general diffusion equation 
in Eq. (I3.1ip . which must reflect, in addition, that within the temporal and spatial reso- 
lution of the variables x and r, the local concentration profile n{r, t) remains uniform and 
stationary. These assumptions can be summarized by the following stochastic equation for 
5n(r + X, t + r) 



dr 



d6n{v + X, t + r) ^^o^^^^ ^^^^ . T f ^^^^^ _ r - r'; r, t] V,, 

(3.18) 

/ (ix2cr""^(| xi - X2 \;t)6n{r + X2,t + r') + f(r + x, t + r). 
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where the function cr~^(| x — x' \;t) is the inverse of the covariance x — x' (in the 
sense that their convolution equals the Dirac delta function), so that its Fourier transform is 
1 / a{k; r, t). The random term f(r + x, t + r) of eq. 03.1 8p is assumed to have zero mean and 
time correlation function given by < f(r + X, t + r)f+ (r + x', t + t') >= 7[x — x', r — r'; r, t], 
with 

7[x-x',r;r,t] = L'Vr,i)Vx ■ y" c/xi6[x - Xi, r; r, t] Vxi<5(xi - x'). (3.19) 
Similarly, the analog of Eq. (12. 6p for the time correlation function C(r, t) is the relaxation 



equation for C(x, r; r, t) = 6n{r + x, t + r)5n(r, t), namely, 
dC{x,T;r,t) _ 0^ 



=D'n{T, t) Vx ■ y^ dr' I rfxi6[x - Xi, r - r'; r, t] Vx, 
J (ix2a'^(xi,X2;t)C(x2,r';r,t). 



(3.20) 



D. Approximate self-consistent closure for the local mobility 6(r, t) 

The generalized theory of non-equilibrium diffusion just presented writes the relaxation of 
the mean value n(r, t), of the covariance cr(r, r'; t), and of the two-time correlation function 
C(x, r;r, t), through Eqs. (13. 6p . (13.80 (or (14.30 ). and (13.200 . in terms of the generalized 
mobility 6[x, r;r,t] or, according to Eq. (I3.15p . in terms of the temporally and spatially 
nonlocal friction function A(^[x, r; r, t]. These equations constitute the general framework 
in which approximations may be introduced to construct a closed system of equations for 
the properties involved. The main purpose of the present subsection is to determine an 
independent closure relation for the local mobility 6(r; t), needed in Eqs. ( 13. 6p and ( 13. Sp (or 
(14. 3p ). in terms of n{r] t)and cr(r, r'; t). 

This, however, will be a relatively involved process. The reason for this is that, according 
to Eq. (I3.17p . the local mobility b{r;t) is an integral of the non-local generalized mobility 
6[x, r;r, t] appearing in Eq. (13.200 for the time-correlation function C(x, r;r, t). Thus, the 
determination of 6[x, r;r, t], is essentially equivalent to the determination of C(x, r;r, t), 
which is intrinsically an involved and rich problem, even under ordinary equilibrium con- 
ditions. Thus, our answer to this problem is equivalent to extending to non-equilibrium 
conditions the equilibrium theoretical approach to calculate these dynamic properties. 

With this aim let us refer to Eq. (I3.2O0 and assume that, within the approximation of 
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local uniformity and isotropy introduced above, C(x, r; r, t) = C(| x |, r; r, t). We then write 
the Fourier transform (FT) of this correlation function as 

C(|x|,r;r,t) = ^ J d'ke~''^'^C{k,T-v,t), (3.21) 

Denoting also the FT of h[\ x |,r;r, t] as b{k,T;r,t), we can then rewrite Eq. fl3.20p in 
Fourier space as 

dC{k,T;r,t) 



dr 



-k^D^n{r, t) dr j b{k, r - r'; r, t)a-\k; r, t)C{k, r'; r, t). (3.22) 



In its turn, the mobility b{k, r; r, t) can be expressed in terms of the FT A({k, r; r, t) of 
AC (I X I, r; r, t) according to Eq. fl3.15p . which in Laplace space reads 

b{k,z;T,t)= 1 + AC{k,z;r,t) ^ (3.23) 

with A^*(A;, z; r, t) = A({k, z; r, t)/C° and with the hat and the argument z meaning Laplace 
transform (LT). Using this result in the Laplace-transformed version of eq. fl3.22p . we finally 
get the following expression for the LT of C{k, r; r, t) in terms of A(*{k, z; r, t) 

C{k, Z; r, t) = k2}on{r,t)X{k;r,t) ■ f^^'^^) 

^ l+AC*(fc,2;r,t) 

Let us notice that we can also introduce the notation by C{k, r; r, t) = n{r, t)F{k, r; r, t), 
with F{k, t; r, t) being the non-equilibrium intermediate scattering function, whose initial 
value F{k,T = 0;r, t) = S{k]r,t) defines the time-evolving spatially local static struc- 
ture factor S{k;r,t). With this more familiar notation it is not difficult to recognize 
in Eq. (13.241) the non-equilibrium extension of the well-known exact expression for the 
LT of the intermediate scattering function in terms of the so-called irreducible memory 
function A(*{k, z;r,t) 1^, 15, 18]. There is, however, a deep fundamental difference be- 
tween this expression for C{k, z;r,t) and its equilibrium counterpart: the initial value 
a{k;r,t) = n{r; t) S {k; r , t) needed in Eq. f l3.24p derives from the nonequilibrium solution 
of the relaxation equation in Eq. (14. 3p . and not from the local equilibrium approximation 
(T'-^ {k;n{r,t)) = [S{k;n{r,t))]^^. Of course, the general expression in Eq. (I3.24p contains 
the conventional equilibrium result as the particular case in which the static structure factor 
S{k]r,t) = a{k;r,t)/n{r;t) is given by its equilibrium value S^'^{k\rf^) = \n'^'^S{k;n'^'^)]~^ . 
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Let us mention that the equihbrium counterpart of Eq. f l3.24p can also be derived without 
appeahng to the phenomenological non-hnear and non-local extension of Pick's diffusion 



equation in Eq. (13.111) . Thus, in Ref. 56| the non-Markovian extension of Onsager's 



theory (referred to there as the "generalized Langevin equation" (GLE) approach) was 
employed to derive the equilibrium version of Eq. (I3.18P , from which the equilibrium version 
of Eq. (I3.24P follows. The value of the phenomenological derivation of the non-linear Pick's 
diffusion equation of Eq. (13.111) is that it is a natural non-linear extension of the more 
rigorously- derived equilibrium linear theory. A similar situation arises when one considers 
the derivation of the result analogous to Eq. (I3.24p for the self component Cs{k, z]r,t) of 
C{k,z;r,t). This result that can also be derived in either of these two manners, both of 
which lead to the following expression for Cs{k, z; r, t) 

Cs{k, z; r, t) = \,^, . (3.25) 

^ ~'~ l+AQ{k,z;r,t) 

In this manner, Eqs. (I3.24p and (I3.25P write the non-equlibrium collective and self time- 
correlation functions C{k, z; r, t) and Cs{k, z; r, t) in terms of the respective irreducible mem- 
ory functions A(*{k, z;T,t) and A(g{k, z;r,t). At this point, with the aim of establishing 
a self-consistent scheme for the calculation of these four properties, we propose to proceed 
along the same lines, and to adopt the same approximations, of the equilibrium SCGLE the- 
ory in its simplest formulation jj]. Thus, we start by adopting the Vineyard approximation 

AC{k,z;r,t) = ACs{k,z-r,t), (3.26) 
along with the factorization approximation 

AC{k, z; r, t) = X{k; r, t) AC{z; r, t), (3.27) 



in which the function X{k; r, t) is a phenomenological "interpolating function" j^, |2l|, given 

by 

\{k;v,t) = (3.28) 

where fcc ~ 27r/ d, where d is some form of distance of closest approach. A simple empirical 
prescription is to choose kc as kc = /cmin, the position of the first minimum (beyond the main 
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peak) of the non-equilibrium static structure factor S{k;r,t) = a{k;r,t)/n{r;t) at position 
r and time t. 

The function A(*{z;r,t) in Eq. 03.271) is the Laplace transform of the r-dependent 
friction function AC*(r;r, t) = AC(r; r, t)/(^0) which can be approximated by the following 
expression 



AC*(r;r,t) = ^ / dk k'' 



S{k;r,t) - 



F{k,T;r,t)Fs{k,T;r,t). (3.29) 



Sik;r,t) 

The derivation of this expression follows, in a first approximation, essentially the same 
arguments employed in the derivation of its equilibrium counterpart, explained in the original 
presentation in Ref. 56| (also reviewed in appendix B of Ref. |2l|). The main aspect 
that needs to be adapted refers to the statistical distribution of the local concentration 
profile of the particles around a particular tracer particle, whose mean and covariance in 
the original derivation refers to the equilibrium distribution, whereas now they refer to the 
mean and covariance of the statistical distribution representing a non-equilibrium state. In 
this manner, the exact results in Eqs. f l3.24p and fl3.25p . complemented with the closure 
relation for the time-dependent friction function in Eq. f l3.29p and the Vineyard and the 
factorization approximations in Eqs. (13. 26 p . (I3.27P and (I3.28p . constitute a closed system of 
equations that must be solved self-consistently. 



IV. FULL NON-EQUILIBRIUM THEORY AND PARTICULAR LIMITS 

In summary, the NE-SCGLE theory is defined in terms of a system of equations for the 
time-evolution of the mean value n(r, t) and of the covariance a{r, r'; t) of the fluctuations 
of the local concentration profile n{r,t) of colloidal particles, namely, 

= D'V ■ 6(r, t)n{T, t) V/3/i[r; n{t)] (4-1) 



and 

da{r, r'; t) 
di 



--D^V-n{r,t) b{r,t)V J dr2£[r,r2;n{t)]a{r2,r';t) 
+ D'V ■ n{r', t) 6(r', t)V j rfrs^ [r', v^- n{t)]a{v2^ r; t) (4-2) 
- 2L)°V ■ n(r, t) b{r, t) V5(r - r'). 
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with £^[r, r';n] = [(5/5/i[r; n]/(5n(r')]. We assume that we can approximate this thermo- 
dynamic matrix as S[r,r + x;rr(t)] ^ S{\x.\;n{r,t)), i.e., by the thermodynamic matrix 
evaluated at a uniform concentration profile with a constant value given by the local and 
instantaneous concentration ?T(r, t) at position r and time t. Then the covariance can also 
be approximated as a{r, r + x; t) ^ cr(|x|; r, t), and the latter equation can also be written 
as 

da{k; r,t) 

(4.3) 



dt 



- 2fc^D°n(r, t)6(r, t)£{k] n(r, t))a{k- r, t) 



where a{k;r,t) = J d^ke'^^'^a{\:x.\;r,t) and S{k;n{r,t)) = (27r)~^ J d^ke~'^^'^S{\x.\;n{r,t)). 

Besides the chemical equation of state (i.e., the functional dependence of the local elec- 
trochemical potential (3fi[r;n{t)] on the concentration profile n{r,t)), the solution of these 
equations require the simultaneous determination of the local mobility function b{r, t) which 
is given, according to Eqs. (13.171) . (I3.23p . and (13.271) . by 



6(r,t) 



1 + 



POO 

/ rfrAC(r;r,t) 

JO 



(4.4) 



The actual calculation of 6(r, t) requires the solution, at each position r and each evolu- 
tion time t, of a system of equations involving the Laplace transform (LT) C{k,z;r,t) = 
drC^k, r; r, t) of C{k, r; r, t) and of its self component Cs{k, r; r, t), as well as the LT of 
the r-dependent friction function A^*(r;r,t), namely, 

cr(/c;r,t) 



C(fc, z; r, t) 



z + 



fc2_D0n(r,t)o-i(fc;r,f) 
1+X{k;r,t) A(*{z;r,t) 



(4.5) 



Csik,z;r,t) 



z + 



k^DO 



(4.6) 



l+A{A:;r,i) AC{z;r,t) 



and 



AC(r;r,t) 



Dn 



3(27r) 



dk k^ 



a{k;r,t)/n{v,t) 



C{k,T;r,t)Csik,T;r,t), (4.7) 



a{k;r,t) 

with X{k;r,t) being the phenomenological "interpolating function" given by Eq. (I3.28p . 

An important aspect of Eqs. (14. 5 1) - (14. 7p above refers to their long-r (or small z) asymp- 
totic stationary solutions, referred to as the non-ergodicity parameters of the corresponding 
dynamic properties. These are given by 

C{k,T;r,t) X{k;r,t)(r{k;r,t) 



/(fc;r,t) = lim 



T^-oo cr(/c;r, t) A(/c; r, t)cr(A;; r, t) + /c2n(r, t)7(r, t) 



(4.8) 
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and 



fs{k; r, t) ^ lim Cs{k, r; r, t) = ^[^'^f,. ,y (4-9) 

where the (spatially and temporally dependent) squared localization length 7(r,t) is the 
solution of 

' - ' ^'^dkk' i"'*:^-'>{!<^-'>:^i'^'f:^^-'' „, (4.10) 



7(r, t) 67r2 [A(fc; r, t)a{k- r, t) + k'^n{r, t)7(r, t)] [A(fc; r, t) + k'^-i(r, t)] ' 

If the solution 7(r,t) of the latter equation is infinite, we can say that at that position r 
and waiting time t the system still remains ergodic, but if it is finite, we say that the system 
became dynamically arrested. 



A. Particular cases and limits 



Eqs. f l4.ip - fl4.10p constitute the full non-equilibrium SCGLE theory. Let us recall that 
C{k,T;r,t) = n{r,t)F{k,T;r,t), with F{k,T;r,t) being the non-equilibrium intermediate 
scattering function, whose initial value F{k,T = 0;r,t) = S{k;r,t) is the time-evolving 
spatially varying static structure factor S{k;r,t) = a{k;r,t)/n{r,t). With this notation, 
Eqs. f l4.ip -f HrT0l) will probably appear more familiar. In fact, it is not difficult to recognize 
in these general equations a number of relevant concepts when adequate limits or cases are 
considered, some of which are discussed in what follows. 

The first obvious general limit to discuss refers to the long evolution-time limit, t — t- oo. 
Assuming static external fields and static thermodynamic constraints, one expects that in 
this limit the solution of Eqs. (14. ip and (14. 3 p will converge to a stationary state, denoted 
by ?T'**(r) and cr'^''(A;; r). This stationary state will be a thermodynamic equilibrium state 
if {dn^'^{r)/dt) and {d(j'^'^{k;r)/dt) vanish due to the fact that the two equilibrium condi- 
tions, V/i[r;n*''] = and S{k;n'^'^{r))o''^^{k]r) = 1, have been attained. Other stationary 
solutions of Eqs. f]4.ip -( l43l) might, however, exist in which the derivatives {dn'^^ (r) / dt) and 
{da^^{k;r)/dt) vanish due to vanishing of the local mobility b{r,t — )■ oo), a condition for 
dynamic arrest. We may, however, disregard the consequences of this second possibility, 
and assume that the system will always be able to reach its thermodynamic equilibrium 
state. Furthermore, let us assume that the system is not subjected to external fields, so that 
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n'*^(r) = fih and a^^{k] r) = nbS{k; /3, n^), with S{k; (3, rib) being the equihbrium static struc- 
ture factor of the homogeneous system. Under these conditions (i.e., full equilibration and 
spatial uniformity), from Eqs. fl4.5p - fl4.7p we recover the equilibrium version of the SCGLE 
theory jl|, Q , and from Eqs. fl4.8p - fl4.10p we recover the corresponding so-called bifurcation 
equations using the terminology of MOT 29 1. 

The full non-equilibrium SCGLE equations, Eqs. fl3.6p -f HrT0|) . can be solved only after 
several elements have been specified. The most basic of them refers to the nature of the 
system, defined by the pair interaction potential u{r), which determines the non-ideal contri- 
bution to the electrochemical potential. This contribution is represented by the term — c[r; n] 
of the chemical equation of state, written as /3/i[r; ri] = +lnri(r) — c[r; n] + /3'?/'(r). The 

functional dependence of c[r; n] on the concentration profile ?T,(r) is a second fundamental 
element that must be specified. One possibility is to propose a theoretical approximation 
for this dependence, which in the language of density functional theory is actually equiva- 
lent to proposing an approximate free energy functional. For example, within the simplest 
approximation, referred to as the Debye-Hiickel or random phase approximation, c[r; n] is 
written as c*-'^^'^-* [r; n] = — /3 j d^v'u{\r — r'|)n(r') (which also defines an approximation for 
the thermodynamic matrix, namely, £'^^^'^\y, r'; n] = 5(r — r')/r2(r) + n(|r — r'|)). 

A third element to specify refers to the external fields and the thermodynamic con- 
straints to which the system is subjected. We have assumed so far that the external fields 
are static and represented by V'(r), whereas the thermodynamic constraints consists of keep- 
ing the temperature field uniform, T{r,t) = T{t) (= l/kBf3{t)), but not necessarily con- 
stant. There is, however, no fundamental reason why we have to restrict ourselves to these 
conditions. In fact, the general equations of the NE-SCGLE theory above can be used, 
within the range of validity of the underlying assumptions, to describe the response of the 
system to prescribed time-dependent external fields V'(r, t) or programmed thermal con- 
straints described by the time-dependent temperature T(t). This would be done by just 
including this possible time-dependence in Eq. (13. 6p through the electrochemical potential 
/i[r, t; n] = fi*(T(t)) + kBT(t) Inn(r) — kBT(t)c[r; n] + ip{r, t). Most commonly, however, we 
assume that such time-dependent fields and constraints could be used to drive the system to 
a prescribed initial state, described by the mean value ^^^(r) and covariance a^{k\ r), for then 
programming the field and the temperature to remain constant afterward, ip{Y,t) = tp{r) 
and T(t) = T for t > 0. The present theory then describes how the system relaxes to its 
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final equilibrium state whose mean profile and covariance are n^'^(r) and a^^{k\ r). 

Describing this response at the level of the mean local concentration profile rr(r, t) is 
precisely the aim of the recently-developed dynamic density functional theory (DDFT) 10, 



60[ • To establish direct contact with this theory, let us consider the limit in which we neglect 



the friction effects embodied in A(^*(r;r,t) by setting b{r,t) = 1 in our main equations, 
namely, Eqs. (14. ip and f l4.3p . We notice that under these conditions Eq. (14. ip corresponds to 
the central equation of DDFT, which has been applied to a variety of systems, including the 
description of the irreversible sedimentation of real and simulated colloidal suspensions 61 1. 



We should also mention that Tokuyama 64, |65| has proposed an equation for the irreversible 



relaxation of n(r, t) which differs from such simplified version of our Eq. (14. ip only in that it 
neglects external forces as well as the effects of the interparticle direct interactions embodied 
in the non-ideal part of the electrochemical potential, i.e., it sets c[r; n{t)] = in Eq. (13. ip . In 
contrast, Tokuyama's theory does include some effects of the direct interparticle interactions, 
as well as of hydrodynamic interactions, on the matrix £[r, r';t] (see Eq. ( 13. 7p ). through 
the replacement of the diffusion coefficient by the short-time self diffusion coefficient 
Ds(n{r;t)) that depends as an ordinary function on the local concentration. Just like 
DDFT, Tokuyama's theory provides a description of the spatially inhomogeneous relaxation 
of the local concentration profile. Furthermore, it seems to predict dynamic arrest for hard- 
sphere dispersions. The current versions of dynamic density functional theory, on the other 
hand, cannot predict dynamic arrest phenomena because of the simplifying approximation 
b{r,t) = 1. 

The theory proposed in the present work shares some elements with both of these the- 
oretical developments, in the sense that it is also aimed at describing the non-equilibrium 
relaxation of the local equilibrium profile. We consider, however, that the description of the 
irreversible relaxation of the macroscopic state of the system is not complete without the 
description of the relaxation of the covariance matrix cr(r,r';t) in Eq. (13. 8p (or (14. 3p ) and 
without the inclusion of the effects embodied in the local mobility function 6(r, t) ^ 1, i.e., in 
the friction function A(*{t; r, t) ^ 0. In this regard it is also important to point out that in 
the limit b{t) — )■ 1 of Eq. (14.30 one can recognize an equation that has been fundamental in 



the description of the early stage of spinodal decomposition [ll|, l66l-l68j . For example, with 
the additional small-wave- vector approximation for S'j-'^{k), namely, Sj'^{k) ^ So+S2k'^+S4k'^, 
this equation is employed in the description of the early stages of spinodal decomposition 
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(see, for example, Eq. (2.11) of Ref. [ll|, in which £4 = 0, or Eq. (23) of Ref. 68]). 

Another particular limiting condition that merits discussion, now in the context of 
the complete theory, corresponds to the quasistatic process, characterized by a trajectory 
rr'-'^-(r,t) and cr''^'(/c; r, t) that satisfies what we refer to as the local equilibrium approx- 
imation. In this idealized process the system is driven from a given initial equilibrium 
state described by rf{r) and cr°(A;;r) to a final equilibrium state described by n'^'^{r) and 
a'^'^{k; r) by extremely slow time-dependent fields and constraints in such a manner that 
{dn'"^' {r , t) / dt) and (9cr''^'(A;; r,t)/dt) virtually vanish due to the fact that at each time t the 
system is allowed to approximately attain the two equilibrium conditions, V/i[r; n^-^-] = and 
a'"'^' {k; r , t) = S^^{k;n'"'^'{r,t)). A quasistatic process, however, is an idealized and rather 
unrealistic concept, at least in the limit of small wave- vectors, in which the relaxation times 
diverge as (see Eq. (14. 3p ). In fact, far more interesting is the opposite limit, in which 
the system, initially at equilibrium with a static field ip^^\r) and temperature T^^\ must ad- 
just itself in response to an instantaneous change of these control parameters to new values 
ijj^-l'\r) and T^-^\ according to the "program" described by V'(r,t) = tP^^\v)9{-t)+ip^f\r)9{t) 
and T{t) = T^^^e{-t) + T^f^e{t) with e{t) being Heavyside's step function. 

Under the conditions described by this instantaneous quench program the predicted non- 
equilibrium trajectory n(r, t) and a{k; r, t) will spontaneously reach the new thermodynamic 
equilibrium state n'^'^{r) and a'^'^{k\ r), unless dynamic arrest conditions arise along this non- 
equlibrium trajectory. This is, of course, the most fascinating possibility, and it was the 
main motivation to carry out the present non-equilibrium extension of the SCGLE theory. 
A simple manner to monitor if this possibility will actually interfere with the process of full 
equilibration is to solve Eq. fl4.10p for the squared localization length 7'^''(r) when n(r, t) and 
a{k]r,t) are given their expected equilibrium values n^'^(r) and cr^'^(A;;r). If the resulting 
value of the dynamic order parameter 7'^''(r) turns out to be finite for r in some portion of 
the system, we should expected the system to become dynamically arrested at least in that 
region. The possible scenarios in which this might be predicted to occur are hidden in the 
full NE-SCGLE equations above and in the specific systems and conditions that might be 
considered. In order to explore how reasonable these expectations may be, in a separate 



work 



12| we apply for the first time the full NE-SCGLE theory above, to the quantitative 



description of the response of a simple model glass-forming colloidal liquid subjected to a 
spatially homogeneous instantaneous quench to conditions where dynamic arrest is expected 
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on the basis of the procedure just outhned. 



V. DISCUSSION AND SUMMARY 

In this paper we have proposed the extension of the self-consistent generahzed Langevin 
equation (SCGLE) theory of colloid dynamics to general non-equilibrium conditions. This 
extension describes in principle the process in which the spontaneous evolution of the system 
towards its thermodynamic equilibrium state could be interrupted by the appearance of 
conditions for dynamic arrest. The main fundamental basis of this general self-consistent 
theory were provided by the general principles of Onsager's theory of equilibrium thermal 



fluctuations, or, better^ 
Markovian conditions 



Dy the extension of Onsager's theory to non- stationary and non- 
57|, whose review was the subject of section |IIl Clearly, this extended 
theory of irreversible processes is in principle applicable to other relaxation phenomena 
outside the realm of colloid dynamics. 

The application of this extension of Onsager's theory to the description of the irreversible 
evolution of the structure and dynamics of a colloidal liquid was carried out in Sect. IIIII The 
resulting non-equilibrium theory of colloid dynamics contains as particular number 
of relevant limiting conditions. For example, the evolution equation for the mean profile 
n(r, t) is found to contain the fundamental equation of dynamic density functional theory 
as a particular limit, whereas the basic equation employed to describe the evolution of the 
static structure factor in the early stages of the process of spinodal decomposition can be 
recognized as a particular limit of the evolution equation for the covariance. The general 
theory, however, also allows its application to the description of the irreversible processes, 
such as aging, associated with dynamic arrest transitions. In particular, it should in principle 
be suitable to describe processes of dynamically arrested spinodal decomposition. 

Let us finally notice that the general non-equilibrium theory of dynamic arrest has built 
in a very natural manner the description of static and dynamic heterogeneities, since at any 
evolution time all the relevant static and dynamic properties are defined at each point in 
space, and cannot a priori be assumed spatially homogeneous. As a zeroth order approxima- 



tion, however, one may simplify the full sel 



-consistent theory assuming spatial homogeneity. 



as it is done in the accompanying paper 121], and for some purposes this simplifying approxi- 



mation may suffice to provide an acceptable first-order scenario of important non-equilibrium 
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processes. 

Since the final value of this general theoretical proposal depends on its actual predictive 
power, in the accompanying paper we illustrate the practical and concrete use of the present 
non- equilibrium theory with a quantitative application to the prediction of the aging pro- 
cesses occurring in a suddenly quenched colloidal liquid, whose static structure factor and 
its van Hove function evolve irreversibly from the initial conditions before the quench to a 
final, dynamically arrested state. As reported there, the comparison of the corresponding 
numerical results with available simulation data seem highly encouraging. 
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Appendix A. FUNDAMENTAL THERMODYNAMIC FRAMEWORK 

This appendix summarizes the essential concepts of the thermodynamic theory of inho- 
mogeneous fluids as a straightforward application of the first and second laws of classical 



thermodynamics 



to a system that cannot be spatially homogeneous. Augmented with 
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58 



62 



63| , the result- 



elementary concepts of the thermodynamic theory of fluctuations 
ing purely phenomenological description involves some basic equations that also appear in 
microscopic statistical mechanical theories, such as density functional theory of inhomoge- 
neous fluids IsQ ]. 

The first law of thermodynamics states that, in a macroscopic system formed by 
particles in a volume V, the total internal energy E is a state function, whereas the second 
law postulates the existence of the entropy S, another state function with the property that 
a closed system will spontaneously search for the state with the maximum S, and this state 
is referred to as the thermodynamic equilibrium state. The functional relationship between 
the entropy and the other extensive variables E, N, and V is referred to as the fundamental 
thermodynamic relation (FTR) of the system, written as 5 = S{E^N ,V). The presence of 
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external fields, however, may cause spatial inhomogeneities in the distribution of matter and 
energy. The description of the possible thermodynamic states of this system then requires 
of more information than that contained in the value of the total properties E, N, and V. 

1. Thermodynamic state space of a non- uniform system. 

For this reason we mentally partition the volume y in a number C of smaller portions 
(or cells), whose internal energy, particle number, and volume, we denote by A^('') and 
V^^\ respectively, with r = 1,2, ...,C. Then, the fundamental thermodynamic relation of 
this system reads S — -SfE, N, V], where E, N and V are C-dimensional vectors with 
components and V^'^^ (r = 1,2,...,C). For the sake of simplicity let us assume 

that the volumes V^^^ arc all equal, V^''~^ = AV = V/C, and remain fixed, so that only the 
variables [E,N] are needed to define a thermodynamic state. Specific values given to each 
of the components N^'^'' of the vector N define a particle number profile, and specific values 
of the components E^'^^ define a specific energy profile E. For notational convenience let us 
also introduce the M-component vector a (with M — 2C) as a = [E, N] . Then, a specific 
particle number profile N and a specific energy profile E define a specific thermodynamic 
profile a. The set of all possible thermodynamic profiles, that we refer to as the entire 
thermodynamic state space T, is then identical to the set of all particle and energy profiles 
that result from giving the components E^^\ N^^^ any value in the range < oo and 

< oo. The fundamental thermodynamic relation of this system, which then reads 

5 = 5[E,N], (1) 

assigns a value of the entropy to any possible thermodynamic profile [E, N] = a e T. 

None of the elements of the state space T of a given system is a priori an equilibrium or 
non-equilibrium state. The second law of thermodynamics can distinguish which element of 
Tis the thermodynamic equilibrium state only after specifying i) the system, ii) the external 
fields acting on its constituent particles and Hi) the global thermodynamic constraints (such 
as isolation or contact with reservoirs) imposed on the system. The system is defined by 
specifying the pair interaction energy u^'^''^"> between two of its particles located at cells r and 
r'. The given array of external fields acting on the particles is described by the corresponding 
total potential energy of one particle at cell r , that we shall denote by ■0^^^. We shall refer to 
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the C-dimensional vector ^ with components -^^^^ (r = 1, 2, C,) as an eocternal potential 
profile. 

The conceptually simplest and most important global thermodynamic constraint that 
may be imposed on the system is total isolation, which prevents the system (of fixed total 
volume V) from exchanging matter and energy with external reservoirs. Thus, the total 
energy E and particle number N are constant, 

c 

E^""^ ^E (= canst.) (2) 

and 

c 

^AtW^TV const.). (3) 

r=l 

We may classify the elements of the entire thermodynamic state space T of a given system 
according to the possible closure conditions, i.e., according to the specific values of {E, N, V). 
Thus, each specific value of {E, N, V) defines a specific subspace t{E, N, V) C T, which 
contains all the thermodynamic profiles a = [E, N] consistent with the referred isolation 
condition. We may then say that any two subspaces t{E, N, V) and t{E', N', V') are disjoint 
unless E = E', N = N', and V = V, and that the union of the subspaces t{E, N, V) for 
all possible values of {E, N, V) is identical to the entire thermodynamic state space T. 
The second law of thermodynamics then states that, for a fixed external potential profile 
^, an isolated system will spontaneously relax from any arbitrary thermodynamic profile 
[E, N] e r{E, N, V) towards the particular profile [Eeq, Ne^] that maximizes the entropy 
within the subspace t{E,N,V). This means that each possible profile ^ will identify a 
member of r(£',iV, as "its" corresponding equilibrium profile [Eeg,NeJ. 

2. Equations of state and conditions for thermodynamic equilibrium. 

The fundamental thermodynamic relation S — <S'[E, N] can also be written in its differ- 
ential form as 

c c 
d5[E,N]/A;B = J]] N](iE('-) - ^/3/x(")[E, NjdArM, (4) 

r=l r=l 

where and denote the functions of the variables E and N defined as 
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For simplicity we may denote by (3 the C-dimensional vector with components f3'^'^\ and by 
the C-dimensional vector with components /J/i*^''-' so that, for example, Eq. (jlj) can also 
be written as rfS'[E, N]/A;b = ^[E, N] ■c/E-/3/i[E, NMN. Eqs. ^ and (ED are, respectively, 
the thermal and the chemical equations of state 58|]. Equilibrium states will satisfy the 
extremum condition dS\Eieq,'^eq\ = which, together with Eqs. ([2]) and ([3]), leads to the 
following set of 2C equations for the 2C variables [Eeg, Neg] 

/3M[E,„N,,] = /3 {= const.) (7) 



/3/i(")[Ee„Ne,] =/3/i {= const) (8) 

for r = 1,2, ...,C. Clearly, these are merely the conditions for internal thermodynamic 
equilibrium, which require that the intensive parameters do not vary from cell to cell. 

3. Thermodynamic theory of fluctuations: covariance and stability matrices. 

The equilibrium value of the thermodynamic profile Req = [Egg, Ngg] is then the solution 
of the extremum condition in Eqs. ([7]) and ([8]). There are, however, instantaneous departures 
from such an equilibrium profile, whose properties can only be described in statistical terms. 
This then means that the thermodynamic profile a = [E, N] must be regarded as a M- 
component random vector, subject to a probability distribution P^''[a] whose mean value 
a is the equilibrium value a'^'^. Thus, we must now recognize that the macroscopic state 
of our system cannot be described simply by indicating the mean value a = ae^; instead, 
it must be described by the full probability distribution function P[a] given, according to 



the thermodynamic theory theory of fluctuations j48|, |58|, |62|, |63| , by the Boltzmann- Planck 
expression P'^'?[a] = exp [(^[a] — 5'[a^'^]) /A;^]- The MxM covariance matrix a'^j = ((5a)(5a)t 
of this distribution function, with elements defined as 

a^J = = P">]iai - af )(a, - af), z,j = 1, 2, M, (9) 

a 

is given by the following exact and general result. 
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eg . geq _ j 



(10) 



where I is the M x M identity matrix and S'^'^ is the equihbrium stabihty matrix, defined as 



At any arbitrary state a (not necessarily an equihbrium state) one can define the second 
differential of theentropy as d'^S[aL]/kB = —da) ■ ^[s]- da, with £^[a] being the M x M matrix 
defined as 



In the last member of this equation, Fi[a\ is the thermodynamically conjugate variable of 
the extensive variable a,, defined as -Fja] = {dS[a\/dai). The conjugate variables -Fj[a] and 
the thermodynamic matrix i^ij[a], are thus defined at any thermodynamic state a. It is, 
however, only when the state a is an equilibrium state that these state functions have an 
important extremum and stability significance. In particular, it is only under conditions 
of thermodynamic equilibrium that the matrix S is the inverse of the covariance matrix, 
according to Eqs. ( ITOl) and ( ITT]) . 

4. Legendre-transformed fundamental thermodynamic relation. 

Just like in ordinary classical thermodynamics, under some circumstances one may prefer 
to express all the previous results not in terms of [E, N] as independent state variables, 
but in terms of [/?, N]. Regarding the internal equilibrium conditions, Eqs. (l7])-(jH]), this 
amounts to eliminate the variables E^'''^ from this set of 2C equations by first solving the 
thermal-equilibrium condition /3^'')[E, N] = (3^ for E, and then substituting the solution, 
denoted as Eig\p^,N], in Eq. ([8]). This leads to C equations for A^W (r = l,2,...,c), 
namely, N] = /3/i(^)[EeJ/3^; N], N] = (r = 1, 2, C), where the functions 

[(3; N] and /^/x'-''"-' [E, N] differ from each other in the set of variables they depend on. This 
procedure is done more formally by defining the Legendre transformation of the fundamental 
thermodynamic relation S = S[E,N], which reads = J'[l3,N] = ^[E, N]/A;b - j3 ■ E . 
The new "thermodynamic potential" N] now plays the role of the entropy S but in 
the thermodynamic state space spanned by the variables [/3,N], and is related with the 
Helmholtz free energy A by J^[(3, N] = —(3 A. 




(11) 




(12) 
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This representation is most convenient under conditions in which the particle system 
is in contact with a thermal reservoir (in our case the supporting solvent) that keeps temper- 
ature constant and uniform, in which the case system is constrained to the thermodynamic 
state subspace T{(3^; N,V) defined by [/3e^,N], with /3e? = (3^ for all cells r. This means 
that within this constrained thermodynamic subspace the FTR can be written in its differ- 
ential form as dj-' = —(3^'/l[(3^; N] ■ dN and the chemical equation of state /i*^'') = n^'^^[/3^; N] 
now expresses the components of the vector /I as a function of the profile N and of the 
thermal reservoir parameter which may then be considered a control parameter. The 
main advantage of this representation is the simplification of the equilibrium conditions in 
Eqs. ([7]), ([H])- Thus, the equilibrium concentration profile N^q is now determined by the 
condition that the electrochemical potential remains spatially uniform, i.e.. 



/3/i(")[Neg;/3] = {= const). 



(13) 



This representation also simplifies the discussion of the fiuctuations (5N around the equi- 
librium concentration profile Ngg. To see this, let us write Eq. ( ITOj) more explicitly as 



(9E 



_ ( 9/3[E,N] 
ON 

( 9^[E,N] 







I 







(eg) 





I 



(14) 



where the subindex "(eg)" means that the thermodynamic derivatives in this equation are 
evaluated at [E, N] = [Eeq,Neq]. By inverting the thermodynamic matrix, along with some 



straightforward thermodynamic algebra, one can show that the covariance 5N5N^ satisfies 



I. 



(15) 



N = Ne 



This result, however, is again Eq. ( !T0|) with a = N and corresponding to the global constraint 
of contact with a thermal reservoir that keeps temperature constant and uniform. 
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